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The University of Missouri-RoUa has 
developed a cloud simulation facility for 
the study of various atmospheric cloud 
processes. The initial relative humidity 
of the air sample put into the cloud 
chamber is a key parameter in virtually 
any experiment and needs to be known 
accurately. This report describes how 
the cloud simulation chamber itself has 
been used as a condensation type hy- 
grometer to calibrate the system's hu- 
midifier. Two distinct and physically 
different methods for inferring mixing 
ratio are used, one exploiting the sensi- 



tivity of aerosol activation to humidity, 
and the other exploiting the sensitivity 
of the rate of growth of cloud droplets 
to humidity. The two methods give 
agreement with each other to within a 
precision of one part per thousand in 
mixing ratio. 
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1. Introduction 

Several countries have developed facilities for 
accurately measuring humidity. Such a national fa- 
cility ideally includes both an ultra stable humidity 
generator, and a primary standard for measuring 
relative humidity. The primary standard is used to 
calibrate the humidity generator, and the humidity 
generator is used to calibrate other portable instru- 
ments which have been sent to the national facility 
for calibration. The gravimetric hygrometer de- 
scribed by Wexler and Hyland [1] utilizes absorp- 
tion of water vapor by a solid desiccant, and 
precision weighing of the absorbed water. The un- 
certainty is about plus or minus one part in a thou- 
sand for the measurement of mixing ratio (mass of 
water vapor per unit mass of dry gas). Gravimetric 
hygrometers are currently the primary standards in 
the United States [2], in Japan [3] and in the United 
Kingdom [4]. A chilled mirror condensation hy- 
grometer [5] is utilized in France as a transfer hu- 



midity standard. This instrument is periodically 
calibrated at the National Bureau of Standards in 
the United States and then returned to France for 
use as a standard in France. The accuracy is about 
plus and minus 0.03 °C in dew point. This corre- 
sponds to an uncertainty in mixing ratio of plus or 
minus 2 parts per thousand, or twice that of the 
gravimetric hygrometer. The chilled mirror hy- 
grometer however, has the great advantage of 
portability. 

The University of Missouri at Rolla has devel- 
oped a cloud simulation facility for studying the 
processes that occur in clouds in the atmosphere. 
Part of the facility is a humidifier, which has an air 
flow rate of 1 L/s, and which produces relative 
humidities close to 100% at temperatures which 
can be set to within 0.01 "C over the range 5 to 
25 °C [6]. Another part of the facility is a cooled 
wall expansion cloud chamber, which produces a 
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supersaturated environment so that water droplets 
can be grown on a sample of cloud condensation 
nuclei (CCN). A third part of the facility is the 
optical system used to observe the Mie scattering 
from the water droplets, thus determining the mean 
size of growing water droplets as a function of 
time. A fourth part of the facility is the aerosol 
generation system, used to supply the CCN sample 
to the simulation chamber. The cloud chamber, op- 
tical system, and aerosol generation system, and 
other components have been described by White et 
al. [6]. 

This report describes how the cloud simulation 
chamber has been used as a condensation type hy- 
grometer to calibrate the humidifier with a preci- 
sion of one part per thousand in mixing ratio. An 
error analysis to determine whether the absolute 
accuracy is as good as the precision has not been 
made. Since an error analysis for the gravimetric 
hygrometer has been made [1], there is more confi- 
dence in the absolute accuracy of the gravimetric 
hygrometer than in our facility. It would therefore 
be desirable to intercompare a gravimetric hy- 
grometer to our facility, but this direct intercom- 
parison wouldn't be feasible since neither facility is 
portable. It would be feasible to compare our hu- 
midity measurement technique to a portable hy- 
grometer such as the chilled mirror instrument 
described by Merigourx and Cretinon [5] and our 
Center would welcome such an opportunity. 



2. Expected Humidity from Humidifier 

The humidifier consists of two columns each 
with 60 vertical glass rods, 0.7 cm diameter, spaced 
on alternate points of a 1.27-cm grid, and enclosed 
in an aluminum cylinder 16.5 cm inside diameter. 
Water flows downward on the glass surfaces, and 
air flows downward in the space between the glass 
rods. The two humidifier columns, each with an 
effective length of 89 cm, are connected in series 
with respect to the air flow. Thus the length of the 
humidifier is in effect 178 cm. The air flow rate is 
1.0 L/s, and the water flow rate is 0.25 L/s. 

The first concern addressed was whether the 
glass rods were isothermal over their length. As 
water evaporates at the rod surface the rods are 
cooled due to the latent heat absorbed. To calcu- 
late the amount of cooling an energy balance is 
performed. To humidify originally dry air the la- 
tent heat flux is 35.6 watts. This heat when supphed 
by the 0.25 L/s liquid water flow produces only a 



0.034 °C temperature drop in the water. Thus the 
glass rods have essentially no axial temperature 
gradients. 

The next concern is how to model the vapor dif- 
fusion in the humidifier. The problem is three di- 
mensional, and would need to be solved 
numerically if one wished an accurate analysis. In- 
stead a crude analysis was performed in which the 
three dimensional problem is approximated as a 
two dimensional problem for which an analytical 
solution exists. Specifically, the humidifier is mod- 
eled as a round tube, having the same hydraulic 
diameter as the cross section of the humidifier. The 
air velocity in the model tube equals the average 
air velocity in the actual humidifier (5.25 cm/s), 
and the length of the model tube equals the total 
length of the two humidifier columns (178 cm). 
The wall of the model tube is wet and isothermal. 
The air enters the model tube dry and at the wall 
temperature. The outlet relative humidity from the 
model tube is sought. 

The solution to the above idealization of the va- 
por diffusion in the humidifier has been given by 
Shah [7]. The hydraulic diameter is 4A/P, where A 
is the cross sectional area of the flow and P is the 
perimeter wetted by the flow. Applied to the hu- 
midifier, A is the area of a 16.5 cm circle minus the 
area of 60 circles of 0.7 cm diameter. P is the cir- 
cumference of 60 circles 0.7 cm diameter. This 
gives a hydraulic diameter of 5.66 cm. The 
Reynolds number based on a 5.66 cm length and an 
air velocity of 5.25 cm/s is about 200, so the air 
flow is laminar. The Peclet number is about 1 15, so 
that vapor diffusion in the flow direction is small 
compared to that in the transverse direction. The 
178 cm length divided by the 5.66 cm length and 
by the Peclet number is a nondimensional quantity 
called jc * by Shah. The mean outlet relative humid- 
ity, called (1 — 0m) in Shah's analysis, is a function 
of only X*. For our x* (0.272) the predicted mean 
outlet relative humidity is 98.45%, equivalent to a 
mean outlet dew point lower than the wall temper- 
ature by 0.256 °C. This result should be treated as a 
crude estimate. The resulting predicted dew point 
deficit of 0.256 "C is an order of magnitude esti- 
mate. It turns out that the observed humidity 
deficit is 0.29 °C, so the agreement is better than 
one could reasonably expect. Shah also gives an 
analytical solution for the equivalent diffusion 
problem between parallel plates. If the humidifier 
is evaluated in terms of equivalent parallel plates, 
the predicted dew point deficit is about an order of 
magnitude smaller (i.e., about 0.03 °C). 
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3. Humidity Measurement in the UMR 
Cloud Simulation Chamber 

The initial relative humidity in the chamber, 
prior to expansion, is a key parameter in virtually 
any experiment. Despite the care that has gone into 
sample humidification, this has proved to be a diffi- 
cult parameter to fix experimentally, and we have 
had to enlist additional means over and above the 
humidifier analysis given above in order to infer 
the humidity's value. 

Two distinct and novel methods of inferring the 
initial mixing ratio from Mie scattering on a 
monodispersed cloud are described below. These 
have proved sufficiently accurate and reliable that 
we have adopted them in all subsequent experi- 
ments. The variable chosen to represent the initial 
water vapor content is mixing ratio, i.e., the num- 
ber of grams of water vapor per gram of dry air 
and is denoted by /-q. 

3.1 Experiment Description 

The experiment begins with a moist aerosol 
laden air sample at temperature T^o and pressure po. 
Usually To is near 20 "C and/?o is near 14.1 psi. The 
aerosol consists of NaCl particles of size near 0.025 
jLtm radius and concentration near 100/cm^ The 
sample's initial relative humidity is typically near 
83%. It is lowered below 100% by immediately 
raising the sample's temperature after it leaves the 
humidifier to avoid condensation loss of vapor (see 
[6]), At a relative humidity of 83%, the aerosol par- 
ticles consist of very small (--'0.05 jLim) solution 
droplets instead of dry particles. The sample is ex- 
panded to give a cooling rate of 10 "C/min, and the 
resulting cloud is observed via Mie scattering of 
laser light to determine drop size as a function of 
time. The aerosol exerts a negligible influence on 
the measurements of interest here. At our smallest 
observable drop size, 0.70 /xm, a 10% change in the 
NaCl particle mass would only change the solution 
drop's equilibrium supersaturation ratio by 
0.000006; an amount negligible in comparison to 
unity. Small variations in the critical supersatura- 
tion of the aerosol w.ould have no appreciable ef- 
fect on droplet growth rates. 

3.2 Multiple Droplet Growth Rate Method 

The condensational growth of water droplets is 
quite sensitive to humidity. Hence droplet growth 
rates can be used to give a measure of the air's 
water vapor content. In our method we take a 



small set of droplet growth rate data containing at 
least two growth rate measurements and covering 
a short period of time. This data set is analyzed 
using drop growth theory to extract two un- 
knowns, mixing ratio and sticking coefficient (fi). 
A whole experiment can be subdivided into numer- 
ous small data sets, and a separate determination of 
ro made for each one. In each case the amount of 
water converted into the liquid state is accounted 
for in the determination of initial (before expan- 
sion) mixing ratio Tq. /-q should of course be con- 
stant. 

The droplet growth rate equation [8] for the i-th 
droplet growth rate measurement can be written 

ro=e/[>A(ri)-H(ai+li)fl/AffPsat.i+'Si*}-'-l]-^ 

-f(4/3)7rA^fli^ (1) 

with e= 0.62 197, p is the total pressure, N is the 
droplet concentration per gram of dry air, e^ is wa- 
ter saturation vapor pressure, T is absolute temper- 
ature, a is the droplet radius, a is the droplet 
growth rate, i)eff is an effective diffusion coefficient 
(see [9]) for water vapor in air, p^^i is the saturation 
water vapor density over a flat surface, 5** is the 
equilibrium supersaturation ratio over a droplet of 
radius a, and 

l=D,i^\aLB/KR,T-{-\p/D), (2) 

\^={l-a/2)K(y-m'jrT/R,y''[ap{y+l)]-\ 



l^=(l//3-l/2)D{2n/R,T) 



1/2 



(3) 
(4) 



Here D is the diffusion coefficient for water vapor 
in air, L is the latent heat of condensation for wa- 
ter, B is the slope of the e^ vs temperature curve, K 
is the thermal conductivity of moist air, Ry is the 
gas constant for water vapor, R^ is the gas constant 
for dry air, a is the thermal accommodation coeffi- 
cient (here we use a= 1), and y is the specific heat 
ratio of moist air. The amount of water in the sys- 
tem is constant, so eq (1) should have the same 
value for all drop growth rate measurement points. 
This fact can be used to eliminate the unknown /8. 
Define a minimization parameter as: 

I i-l 



2 .2 (roi-/-oj)V[8(roi-roj)]^ 



x^=^-^ 



2 [8(roi-roj)]- 

ij 



(5) 



where 8(roi— roj) denotes the uncertainty in the 
knowledge of roi— ^oj- This uncertainty is calculated 
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based on a 0.05 jam uncertainty in a, a 10% uncer- 
tainty in a, and a 5% uncertainty in A^. The con- 
densation coefficient enters into x^ through 1. Since 
ro is constant x^ should be zero or at least small. 
Using the computer routine STEPIT, ji can be 
varied to minimize x^- This yields a value for ^ 
which can be put back into eq (1) to get ro. This 
process can be repeated for each data subset taken 
during the experiment to give multiple measure- 
ments of the same quantity Tq. 

3.3 Cloud Arrival Time Method 

During a given expansion, the cloud is observ- 
able soon after saturation (100% relative humidity) 
is reached, and this observation provides a good 
determination of the saturation event. The cloud 
droplets can be observed starting at 0.70 fxm radius. 
The small time increment, on the order of one sec- 
ond, between saturation and first cloud observa- 
tion, can be accounted for by droplet growth 
modeling. For the special case of expansions using 
constant cooling rates (temperature is a linear func- 
tion of time), this time increment is very insensitive 
to /"o. The method simply involves observing the 
cloud arrival time, subtracting the above time in- 
crement from this to determine the time at which 
the gas sample reached 100% relative humidity, 
noting the gas temperature and pressure at this time 
(100% relative humidity), and finally calculating 
the mixing ratio ro of the gas from this information. 

This method is only applicable to experiments 
using constant cooling rates which make the above 
time increment independent of ro. Also the method 
is applied only to experiments using monodispersed 
CCN aerosols and fairly fast expansion rates. These 
conditions lead to monodisperse clouds (all drops 
are the same size), and allow the use of Mie scatter- 
ing of laser light as the drop sizing method [6]. 

The time increment between the 100% relative 
humidity point and the cloud first observation 
point is calculated from droplet growth theory, 
and does depend on the sticking coefficient. How- 
ever this particular situation involves freshly pro- 
duced droplets, and fresh water surfaces are 
thought to have large sticking coefficients near 
unity [10]. The sticking coefficient results taken 
from the droplet growth rate method described 
above confirm this for the early stage of the cloud's 
lifetime. Hence, in the time increment calculation, 
a sticking coefficient of unity is assumed. The influ- 
ence of Sticking coefficient is modest here in any 
case. A change in sticking coefficient by a factor of 



2 (or 5) leads to a change in ro by an amount 
0.000006 (or 0.000023). 

The cloud arrival time method hasn't been opti- 
mized and therefore improvements could be 
achieved. For instance, the influence of the sticking 
coefficient could be minimized by choosing a 
slower expansion rate and a more sharply 
monodispersed CCN aerosol. This would reduce 
the change in temperature, pressure, and humidity 
during the time interval during which the cloud 
drops grew to observable size. 

3.4 Numerical Cloud Model Tests 

As discussed above, a numerical cloud model, 
which calculates gas thermodynamics and water 
droplet growth processes, is used in this analysis. 
For this we use the numerical model presented by 
Hagen [9]. This model has been successfully tested 
against older results in the literature [9] and against 
the NASA Analytic Simulator [11,12]. Anderson, 
Hallett, and Beesley [13] presented an extended so- 
lution to the droplet growth problem and gave an 
intercomparison among the leading droplet theo- 
ries. They conclude that Carsten's [8] solution of 
the droplet growth problem is quite adequate. This 
is the solution on which our cloud model is based. 
Furthermore Anderson et al. [13] presented numer- 
ical results from their extended model for numer- 
ous test cases. We calculated these cases with our 
model and agreed with their results to within 2% 
for drop size. 



4. Results 

The experimental results from five different days 
of work are given in table 1. The run number gives 
a unique identification number for each experi- 
ment. The first six digits of the run number give the 
date (month-day-year) on which the experiment 
was performed. These experiments run from Sep- 
tember 1986 to November 1987. The column la- 
beled ro(MDG) gives the values for the mixing 
ratio obtained from the multiple drop growth rate 
method. This represents the average of repeated 
measurements of this quantity during the course of 
the experiment. Since ro represents the initial value 
of the mixing ratio, it is fixed and all of these re- 
peated measurements should yield the same value. 
cr(MGR) denotes the standard deviation of the in- 
dividual results around the average and is a mea- 
sure of experiment consistency. Note that cr(MGR) 
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Table 1. Numerical cloud model tests 



Run 


/•o(MGR) 


o-(MGR) 


ro(CAT) 


Aro 


091086.04 


0.012582 


0.000008 


0.012591 


0.000009 


091086.05 


0.012574 


0.000012 


0.012550 


-0.000024 


091086.06 


0.012580 


0.000022 


0.012570 


-0.000010 


091086.07 


0.012605 


0.000005 


0.012615 


0.000010 


111886.A1 


0.012618 


0.000010 


0.012630 


0.000012 


111886.A3 


0.012733 


0.000007 


0.012752 


0.000019 


111886.A4 


0.012551 


0.000010 


0.012565 


0.000014 


121786.01 


0.012480 


0.000005 


0.012488 


0.000008 


121786.02 


0.012483 


0.000007 


0.012485 


0.000002 


121786.03 


0.012503 


0.000006 


0.012512 


0.000009 


121786.04 


0.012521 


0.000013 


0.012533 


0.000012 


121786.05 


0.012509 


0.000006 


0.012519 


0.000010 


121786.06 


0.012520 


0.000008 


0.012524 


0.000004 


121786.07 


0.012530 


0.000011 


0.012526 


-0.000004 


121786.08 


0.012523 


0.000008 


0.012541 


0.000018 


121786.09 


0.012525 


0.000012 


0.012523 


-0.000002 


121786.10 


0.012528 


0.000045 


0.012524 


-0.000004 


080787.01 


0.012668 


0.000018 


0.012689 


0.000021 


100687.01 


0.012510 


0.000008 


0.012506 


-0.000004 



is usually small, indicating good consistency be- 
tween the various measurements taken in a given 
run. /"o(CAT) denotes the initial mixing ratio as de- 
termined by the cloud arrival time method. Aro de- 
notes the difference between the two methods, i.e., 
Aro=ro(CAT)-/-o(MGR). 

The average of the absolute difference between 
ro(MDG) and ro(CAT) is 0.000010, with standard 
deviation around the average of 0.000007. This av- 
erage difference between the two methods corre- 
sponds to a difference of 0.012 "C in dew point. 
Hence the two methods are in good agreement 
with each other. The average difference between 
measured humidifier temperature and dew point is 
0.29 "C, corresponding to a relative humidity of 
98.1% for the air coming from the humidifier. This 
figure is in agreement with the rough estimate of 
expected humidifier performance under these con- 
ditions (0.26 °C difference between humidifier tem- 
perature and dew point, or 98.45% relative 
humidity for air leaving the humidifier) that was 
given above. 



5. Conclusion 

The problem of determining the water vapor 
content of gas sample is addressed using the con- 
densational growth of airborne water droplets as 
the measurement tool. Two distinct methods are 
employed. One is based on multiple measurements 
of droplet growth rate; the other is based on the 



fact that cloud is observable almost immediately 
after the gas sample is brought through 100% rela- 
tive humidity. These two quite different methods 
give good agreement with each other. The average 
difference in calculated initial mixing ratio is 
0.000010 grams of water per gram of dry air (corre- 
sponding to about 8 parts in 10*), with a standard 
deviation of 0.000007 gram/gram-air. This corre- 
sponds to a dew point difference of 0.012 °C. This 
agreement gives us substantial confidence in our 
determination of water vapor content. 
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